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Abstract 

We present a quantitative description of J /tp and ip' suppression in proton- 
nucleus and nucleus- nucleus collisions at CERN energies. We use a conventional 
hadronic framework based on nuclear absorption plus final state interaction of the 
J/tj) or tp' with co- moving hadrons. 
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Charmonium suppression due to Debye screening in a deconfined medium was pro- 
posed in 1986 by Matsui and Satz and found experimentally by the NA38 Collabo- 
ration 0. However, it was claimed very soon |3], [|, |5|] that this phenomenon, which is 
also present in pA collisions, could be due to the absorption of the pre-resonant cc pair 
in the colliding nuclei. It is nowadays known ^ 0] that the J/ip data on pA and on AB 
collisions with a light projectile can indeed be described by nuclear absorption with an 
absorptive cross-section (Jab s — 7.3 ± 0.6 mb. 

Recently the NA50 Collaboration has found an anomalous J/ip suppression in PbPb 
collisions, i.e. a suppression which is substantially stronger than the one obtained from 
nuclear absorption with the above value of a a b s [0, || • Two different interpretations of this 
anomalous suppression have been proposed in the literature. One is a scenario in which 
there is an extra suppression of the J/ip due to interactions with co-moving hadrons 
(co- movers) ||, [II], |TTJ. The other interpretation [|TJ, [13], [L4j assumes that when the local 



energy density is larger than some critical value (taken to be around the maximal one 
reached in a central SU collision), there is a discontinuity in the J/ip survival probability. 
(See also \T5\ for ideas based on percolation of strings.) 



In recent papers Kharzeev, Lourengo, Nardi and Satz || and Vogt |16| have claimed 
that a quantitative analysis of the data allows to conclude that the co-mover model 
cannot describe the data, whereas a quark-gluon plasma interpretation describes them 
well. In the present work we re-examine this point by reanalyzing all available data, 
using the final 1995 NA50 results ||, in the conventional co-mover scenario mentioned 



above (see also KIT 



Nuclear absorption: We describe it in the probabilistic model of Ref. [[|. Let us con- 
sider first proton-nucleus collisions. In this model, the pre-resonant cc pair is produced 
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at some point z inside the nucleus and scatters with nucleons on its path at z' > z, with 
an absorptive cross-section a abs . This produces a change in the A dependence of the J /if) 
inclusive cross-section. For nucleus-nucleus collisions this change, at impact parameter 
b and transverse position s, is given || by 

Q abs (h s _ [1 - exp(-A T A (s) a abs )}[l - exp(-E T B (b - s) a abs )} 

{ , ) ~ <?L ABT A (s)T B (b-s) • (1) 

Here T A and Tg are the nuclear profile functions normalized to unity, determined from 
a standard Saxon- Woods density p(r) = Po/(l + ex P [( r — Ra)I°\)i w ith R A = 1.14 A 1 / 3 
fm and a = 0.545 fm []n|. a abs is the absorptive cross-section. In the following we take 
&abs = 7.3 ± 0.6 mb which gives the best fit to the pA data [BJ. Note that in Refs. 0. E| 
a smaller value, a abs = 6.2 ±0.7 mb, has been obtained. This is due to the fact that the 
approximate expression of nuclear absorption in Ref. || was used, instead of Eq. (H). 
Also note that S abs = 1 for a abs = 0, so expression ([!]) has the meaning of a survival 
probability of the J /if) due to nuclear absorption. 

Since we are aiming at a quantitative analysis it should be emphasized that the 
probabilistic formula, with its longitudinal ordering in z, can only be true in the low 
energy limit. Therefore it is important to evaluate the uncertainty resulting from using 
this formula at yfs ~ 20 GeV. In a recent paper [HJ the equivalent of Eq. ([!]) has been 



derived in a field theoretical approach. The obtained formula is valid at all energies and 
coincides exactly with ([!]) in the low energy limit. It is amazing that the differences 
between the results obtained with this exact formula and the ones obtained from Eq. (|TJ) 
are less than 1 %. 

Absorption by co-moving hadrons: The survival probability of the J /if) due to ab- 
sorption with co-moving hadrons is given by (see || [|, |K| and references therein) 
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S co (b, s) = exp 



This formula is obtained assuming longitudinal boost invariance of hadronic densities 
and isoentropic longitudinal expansion (i.e. a decrease of densities with proper time in 
1/t). Transverse expansion is neglected. Ny°(b,s) is the density of hadrons per unit 
transverse area d 2 s and per unit rapidity at impact parameter b. All species of hadrons 
are included in Ny°. (If we consider only the process i[) + p—tD + D + -- - as in Ref. 
0, the value of Ny°(b,s) has to be decreased and that of a co increased by the same 
percentage amount.) 

In order to have a smooth onset of the co- movers and to avoid any threshold effect, 
it is natural to take for Nf the density of hadrons per unit rapidity in a pp collision, i.e. 
Nf = [3/(irRp)] dN~ /dy(y* = 0) ~ 1.15 fm -2 . This coincides with the value introduced 
in Ref. 0. Because of this choice of Nf, the ^-function in Eq. (0) is irrelevant; we 
hace checked it numerically. Thus, Nf cannot be regarded just as a free parameter. 
Moreover, small changes in the value of Nf can be compensated by smaller changes in 
cr co , without spoiling the quantitative comparison to the data. The argument of the 



log is the interaction time of the J/ip with co- moving hadrons. In Ref. [10] a different 
expression for the interaction time based on interferometry radii was used. The results 
obtained with the two expressions are practically identical. More precisely the expression 
(7) in Ref. fllPf , with the initial time To = 1 fm/c used there, gives practically the same 
results as our Eq. (|2|) with Nf = 1.15 fm~ 2 . 

a co is the co-mover cross-section properly averaged over the momenta of the colliding 
particles (the relative velocity of the latter is included in its definition) and over the dif- 
ferent species of secondaries. Unfortunately, the value of a co is not known experimentally. 
This is, of course, the main limitation of the co-mover scenario. Different theoretical cal- 



dilations of the J/^-hadron cross-section based on the multipole expansion in QCD p0| 
differ from those which include other non-perturbative effects jn|] by at least a factor 20 
for -y/i ~ 5 GeV. Other references [2^] obtain values for the J / ip-N cross-section at high 
energy of 4 6 mb and a ratio cr^'~ N / 'a J l^~ N ~ 3 -f- 4 in agreement with geometrical 
considerations. Our value crf^/af Q = 10 (see below) is much larger than its asymptotic 
value. This is consistent with the very different behaviour of the two cross-sections near 
threshold. Note, however, that Eq. (§) is the result of an integration from time r to 
freeze-out. For times close to r , one is dealing with a dense interacting parton system 
and thus the precise relation between a co and the J/?/>(?/>')-hadron cross-section is not 
established. In this situation, inverse kinematic experiments could help to determine 
the actual role of co-movers in J /if) suppression. Phenomenologically, the value of a co 



obtained here allows to make predictions at other energies, in particular for RHIC p3 
Note that S co (b, s) = 1 for o co = 0. The effects of the co-movers in proton-nucleus 

collisions turn out to be negligibly small. 

The inclusive cross-section for J/ip production in nuclear collisions is then given by 



4 B (b) = — [d 2 s m(6, 8 ) S abs (b, s) S co (b, s) , (3) 

where 



m(b, s) = AB a pp T A (s) T B (b - s) . (4) 
We use a pp = 30 mb. From Eqs. (|j) and ([|) we see that for Drell-Yan pair production 

{(Tabs = <3co = 0), I®B = AB I^jn- 

In the dual parton model (DPM), Ny°(b, s) is given by |24], |25[ 



N v°( b , s) = [Nx m A {b, s) + N 2 m B {b, b-s) + N 3 m(b, s)] 9{m B {b, b - s) - m A (b, s)) 
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+ [N[ m A {b, s) + N' 2 ttib (6, b - s) + N' z m(b, s)]9 (m A (b, s) - m B (b, b - s)) 



(5) 



Here m is given by Eq. @) and m A , tub are the well known geometric factors [25, |2~6 l 



m A (B)(b, s) = A(B) T A{ b)(s) 1 - exp {-a pp B(A) T B{A) (b - s 



(6) 



The coefficients iVj and N- are obtained in DPM by convoluting momentum distribution 



functions and fragmentation functions |24]. Their values (per unit rapidity) for the 
rapidity windows and energies of the NA38 and NA50 experiments are given in Table 1. 
The rapidity density of hadrons is given by 



dN c 



— [d 2 b { d 2 s N™(b,s) 

r AR J J 



(7) 



dy o AB 

with a AB = fd 2 b (1 - exp [—cr pp ABT A B (b)]), T AB {b) = J d 2 s T A (s)T B (b - s). 

The obtained densities of negative hadrons at y* = for pp, SS, SAu and PbPb are 
compared in Table 2 with available data |27], |28], ^9fl , using in each case the centrality 
criteria (in percentage of total events) given by the experimentalists. At this point we 
would like to comment on the differences between DPM and the scaling in the number 
of participants known as wounded nucleon model (WNM, for a review see e.g. f30fl). 
The former gives a multiplicity that increases faster with centrality. This is due to the 
presence of strings of type quark-antiquark |23| . Since these strings contribute only at 
midrapidity, the difference between the two models is maximal at y* ~ 0, where the 
NA38/NA50 spectrometer is located, and quite small at the negative values of y* of the 
NA50 E T calorimeter. The data of the NA35 @ and NA49 Collaborations on the 
rapidity distribution of negatives in central SA and PbPb interactions show an agreement 
with the WNM in the fragmentation region and an excess of 20 4- 30 % at y* ~ 0, as 



expected from DPM. We have checked 123 ] that the correlation between Ex and Ezdc 
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(the energy in the NA50 zero degree calorimeter) in PbPb has a small concavity in DPM, 
resulting in a fit to the measured correlation which is at least as good as the straight line 
obtained in the WNM. 



E T dependence: To determine the Et dependence we need the Et — b correlation, i.e. 
the Et distribution at each impact parameter P(Et, b). The ^-dependence of the ratio 
J/i/j over Drell-Yan is given by 

mp ,_ fd%P(E T ,b)I J A ^(b) 

where I^ B (b) is given by Eq. (Bf) and I^g (b) is obtained (up to a normalization constant) 
from Eq. (BJ) with a abs = cr co = 0. 

It is clear from the discussion below Eq. (I7|) that, in the region of the NA50 calorime- 
ter, P(Et, b) in DPM is very similar to the WNM one || j31fl . However, for consistency, 
we are going to use the DPM distribution: 



P(E T , b) = . 1 exp 
l2irq 2 aN™(b) 



[E T -qN c °{b)Y 



(9) 



2q 2 aN™(b) 

where q and a are free parameters and Ny°(b) = J d 2 s N^°(b, s) is obtained from @ with 
the coeficients iVj and N[ obtained in DPM. They correspond to the density of neutral 
particles in the rapidity windows of the NA38 (SU) and NA50 (PbPb) calorimeters. 
The parameters q and a are obtained from a fit of the Et distributions for dimuon pair 
production above the J /if) mass. The resulting fits are quite good |23| and give: q = 0.65 
GeV, a = 1.5 for SU and q = 0.78 GeV, a = 1.5 for PbPb. It is interesting to note 
that the value of q for PbPb is identical from the one obtained from the best fit to the 
Et — Ezdc correlation. The value of a is poorly determined but affects very little the 
results below. A value a = 1 is also consistent with the data. Our value a = 1.5 is 
the one expected from a Poissonian distribution of clusters (resonances) with an average 



cluster multiplicity of 1.5. This value agrees with the one obtained when clusters are 
identified with a realistic mixture of known resonances and direct particles |32| . 



Numerical results: We present the results for J/ip and ip' suppression obtained with 
two sets of parameters. Set I corresponds to nuclear absorption alone: cr a b s = 7.3 mb and 
°to = a to = 0- Set II contains the effect of the co-movers: a a b s = 6.7 mb, af a = 0.6 mb 
and afo = 6.0 mb (Nj = 1.15 fm -2 as discussed previously). The absolute normalization 
(corresponding to a^ p in (J/ip)/AB and to a^ 1 jo^Y in ip^/DY, in the acceptance of 
the NA38 and NA50 experiments), is a free parameter which, for each Set, has been 
determined from a best fit to the data. 

The results for J j ip suppression versus AB are presented in Fig. 1. Nuclear absorption 
alone, Set I, gives a \ 2 /dof = 1.1; although this value is quite good, the experimental 
PbPb point lies well below the theoretical curve for this Set (by ~ 3 standard deviations). 
However, Set II gives a satisfactory description (x 2 /dof = 0.2) of all points. We also see 
that the effect of co-movers is much smaller in SU than in PbPb. 

We turn next to the Et dependence. Using Eqs. @ and @ we compute the ratio 
R(Et) for SU and PbPb in the five Et intervals of the NA38 and NA50 experiments. 
In order to exhibit all the results in the same figure we plot the ratio R versus L. 
This variable is a measure of the centrality of the collision. The average value of L in 
each Et bin is given in the experimental papers ||. However, the value of L is largely 
irrelevant since we are comparing the measured suppression in specific Et bins with 
the model calculations in the same Et bins; it only provides a scale for the horizontal 
axes. For consistency, we have to take, for each Et, the same value of L used by 
the experimentalists. Note that the first calculations of L by NA50 [F[ used a sharp- 
surface approximation for the nuclear density. More recent calculations are based on 
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a Saxon- Woods density and are in better agreement with other calculations available in 
the literature (e.g. p3| ). 

The results for Sets I and II for J/ip suppression ((J/ip)/DY) are given in Fig. 2 for 
all pA and AB data as a function of L. In Fig. 3 the same results for SU and PbPb are 
presented as a function of Et in the form of a continuous line (not as the average in each 
Et bin as in Figs. 2 and 4). Set I gives a x 2 /dof = 27.9 for all pA, SU and PbPb data, 
indicating that nuclear absorption alone fails very badly. On the contrary, without the 
PbPb data the best fit with Set I gives x 2 /dof = 0.9. Set II gives x 2 /dof = 2.7 with 
only pA and SU data. So without the PbPb data it is hardly possible to decide whether 
co-movers are present or not - although the x 2 /dof is better without co-movers. When 
PbPb data are included Set II gives x 2 /dof = 4.3. What prevents this value of being 
smaller is the peculiar L or Et shape of the PbPb data which cannot be reproduced in 
our simple approach. 

The results for ip' suppression (tp'/DY) are presented in Fig. 4. Set I gives x 2 /dof = 
14.3 and Set II x 2 /dof = 1.3 for all pA, SU and PbPb data. 

The treatment of the co-movers presented above is similar to the one in Ref. ||. 
However, it differs from it and from previous treatments §|, |10| in that we use the DPM 
expression for the density of hadrons, Eq. ([5]), instead of assuming it to be proportional 
to either the number of participants or to Et- Another difference resides in the nuclear 



densities. In the present work (and also in |T2[) calculations have been done with the 
nuclear density described after Eq. (|l]), whereas in || [16| the 3-parameter Fermi distri- 
bution of Ref. |31| is used. Using the latter and keeping all other parameters as above, 
we have obtained a J/ip suppression between the first and the last Et bins which is 7 % 
larger in SU and 4 % larger in PbPb. 

In conclusion, the data on both J/ip and ip' suppression can be described in a co- 
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mover approach with a small number of free parameters, which take reasonable values. 
In this approach there is no discontinuity in any observable. One obtains a monotonic 
decrease of the J/ip and ip' over Drell-Yan ratios from the most peripheral to the most 
central collisions. A clear departure from such a behaviour would rule out the co-mover 
description of J/ip and ip' suppression presented above. It is also important to compute 
the J / ip suppression at RHIC in the two approaches using the parameters determined at 
SPS. This suppression is expected to be quite large and will possibly be wildly different 
in the two scenarios. 
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Table captions: 



Table 1. Coefficients (per unit rapidity) in Eq. ([5]) in the rapidity windows of the 
electromagnetic calorimeters (corresponding to neutral multiplicity) and dimuon detector 
(corresponding to multiplicity of charged plus neutrals) and at energies of the NA38 and 
NA50 experiments. When using these values in Eq. (Q) one should put A = A pro j ect u e 
and B = A target . 

Table 2. Negative particle densities at midrapidity obtained with Eqs. (^) and (0) (th), 
compared to experimental data |27|, |28|, (exp) . Percentages of total events (given in 
the experimental papers), and corresponding impact parameters considered, are shown. 



14 



Figure captions: 



Figure 1. J/ip suppression versus AB: Set I (dotted line) and Set II (solid line) com- 
pared to the experimental data ||. The normalization factors are 2.01 nb/nucleon 2 for 
Set I and 2.08 nb/nucleon 2 for Set II. Note that the calculations have been performed 
only for those nuclei where data exist. The obtained values have been joined by straight 
lines. 

Figure 2. Ratio (J/ijj)/DY (Eq. (§)) versus L (fm): Set I (dotted line) and Set II (solid 
line) compared to the experimental data M. Results obtained as an average over each 
experimental Et bin have been joined by straight lines. The normalization factors for 
the theoretical lines (giving the \ 2 /dof indicated in the text for all pA, SU and PbPb 
data included in the fit) are 38.32 for Set I and 45.50 for Set II. 

Figure 3. Ratio (J/ip)/DY (Eq. (||)) versus E T (GeV) compared to the experimental 
data P| , in the form of a continuous line (not as an average over each Et bin as in Figs. 2 
and 4), for SU (upper figure) and PbPb (lower figure). Conventions and normalizations 
are the same as in Fig. 2. 

Figure 4. Ratio tp'/DY (Eq. (g)) versus L (fm): Set I (dotted line) and Set II (solid 
line) compared to the experimental data M. Results obtained as an average over each 
experimental Et bin have been joined by straight lines. The normalization factors for 
the theoretical lines (giving the y 2 /dof indicated in the text for all pA, SU and PbPb 
data included in the fit) are 0.299 for Set I and 0.723 for Set II. 
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Table 1 







N 2 


N 3 


K 


N> 2 




SU (-1.2 < y* < 1.2) 


0.2096 


0.2746 


0.1598 


0.2827 


0.2015 


0.1598 


P6P6 (-1.8 < y* < -0.6) 


0.3549 


0.0548 


0.0946 


0.3198 


0.0899 


0.0946 


5C7 (0.0 < y* < 1.0) 


0.8433 


0.6003 


0.4995 


1.0854 


0.3582 


0.4995 


PbPb (0.0 < < 1.0) 


0.5891 


0.8086 


0.4248 


0.3685 


1.0292 


0.4248 



Table 2 



Reaction 




dN-/dy\Z£* 


pp 


0.73 


0.76 ±0.04 


SS (11 %, b < 2.7 fm) 


19.3 


19.0 ± 1.5 


SAu (1.3 %, b < 1.3 fm) 


56.3 


59.0 ±3.0 


(5 %, b < 3.4 fm) 


207 


195 ± 15 
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Figure 1 
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Figure 2 
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Figure 3 




Figure 4 
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